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Abstract 

We pursue the analysis of a set of generalized DGMLY sum rules for the 
electromagnetic chiral parameters at order e 2 p 2 and discuss implications for 
effective Lagrangians with resonances. We exploit a formalism in which charge 
spurions are introduced and treated as sources. We show that no inconsistency 
arises from anomalies up to quadratic order in the spurions. We focus on the 
sum rules associated with QCD 4-point correlators which were not analyzed in 
detail before. Convergence properties of the sum rules are deduced from a gen- 
eral analysis of the form of the counterterms in the presence of electromagnetic 
spurions. Following the approach in which vector and axial- vector resonances 
are described with antisymmetric tensor fields and have a chiral order, we show 
that the convergence constraints are violated at chiral order four and can be 
satisfied by introducing a set of terms of order six. The relevant couplings get 
completely and uniquely determined from a set of generalized Weinberg sum- 
rule relations. An update on the corrections to Dashen's low-energy theorem 
is given. 
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1 Motivation 



Virtual photons play an important role in low energy QCD for isospin breaking phenom- 
ena and must also be taken into account in connection with precision experiments. The 
formalism for extending the framework of the standard ChPT^ |5J to include virtual 
photons at one-loop was developed, some time ago, by Urechjlj. Several extensions were 
later performed to accommodate leptonsfSJ, to the weak non-leptonic sector [Hj, to the 
anomalous sector jjj as well as to the baryon sector [Ej- The problem of the determina- 
tion of the new set of low-energy coupling constants (LEC's) which appear in Urech's 
Lagrangian is of obvious practical importance in relation with radiative correction cal- 
culations. Many such calculations have been performed recently: low energy tttt or ttK 
scattering EEE EU El El EH ESI pionic as well as pi-kaonic atoms (El E3 EH1 , the set 
of semi-leptonic decays of the KaonfTHl 120] . processes from the anomalous sector (2*1*1 I7j. 
the problem of computing the hadronic contribution of the muon g — 2 from r decay 
data [22]. the theoretical prediction of the CP-violation ratio e'/e for the Kaon [""3*1 I24| I25j 
to mention just a few illustrative examples. Another important application of the elec- 
tromagnetic low-energy couplings (LEC ) is to the question of the chiral corrections to 
Dashen's low-energy theorem |26j and the determination of the light quark masses (e.g. 
|27j). Interesting questions were raised recently concerning the proper definition of the 
quark masses in the presence of QED|28|. 

The problem of the determination of the EM LEC's was addressed in several papers 
some time ago |29l I5UI I31j but the results remain incomplete and sometimes contradictory. 
In this paper, we continue the discussion started in ref. [""""" (hereafter referred to as (I)). 
In (I) it was shown that the electromagnetic LEC's obey integral sum rule representa- 
tions which generalize the classic DGMLY sum rule [""2] f° r the 7r + — 7r° mass difference. 
The integral representations have the form of convolutions involving a pure QCD n-point 
Green's function with n = 2, 3 and 4 together with the free photon propagator. These 
representations serve a number of purposes. For instance, they can be used to recover 
the chiral scale dependence of the LEC's, to determine how these depend on the £ gauge 
parameter and, also, by studying the convergence properties, to determine which LEC's 
are affected by short distance ambiguities. Finally, these representations can be used as a 
method to provide approximate determinations of these couplings. 

This was exploited in (I) to discuss the set of couplings which are associated with 2- and 
3-point Green's functions. The method consists in constructing rational approximants to 
the Green's functions, which is in the spirit of the large N c expansion |33). The parameters 
of the rational functions are constrained by the chiral Ward identities, by physical input 
on resonance masses and couplings and by short-distance matching conditions. In the 
present paper, we address the more complicated case of the 4-point functions. The way of 
insuring that the Ward identities are satisfied is to start from a chiral Lagrangian including 
resonances. Such a Lagrangian was constructed in ref."""I] which uses the antisymmetric 
tensor formalism |2j to describe the vector and axial- vector resonances. This formalism is 
particularly convenient for discussing the ordinary LEC's because one can assign a chiral 
order to each resonance field such that all the terms which are of order four, which include 
one multiplet of vector, axial-vector, scalar and pseudo-scalar resonances were listed in 
ref. [HI]. It was shown in a subsequent paper (HE] that if one imposes a number of short 
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distance matching conditions, the determination of the 0(p 4 ) LEC's becomes independent 
of the specific representation of the resonances in the chiral Lagrangian. 

The Lagrangian of ref.[31] was employed by Baur and Urechj2H] (BU) to discuss the 
electromagnetic LEC's. For this application, however, one realizes by using the integral 
representation that the Green's functions are needed under kinematical conditions which 
are different, in general, from those considered in[SS] such that the Lagrangian of ref. |34| 
is no longer sufficient to ensure the proper short distance constraints. This deficiency 
manifests itself in the fact that the results of BU for the electromagnetic LEC's fail to 
satisfy the correct scale dependence, which is an effect of leading order in N c . In this 
paper, we will show that a simple extension of the Lagrangian |34j to include a set of terms 
which are of chiral order six allows one to satisfy the short distance constraints in the chiral 
limit which are necessary for a consistent determination of the electromagnetic LEC's and 
which are associated with the 2- and 3-point functions considered in [SJ as well as the 
4-point functions introduced in this paper. Finally, we will present estimates for the six 
LEC's which are associated with these 4-point functions. 

2 Chiral calculation of the basic correlators 
2.1 Definition of the basic correlation functions 

In the presence of a dynamical photon field it is convenient to extend the source part 
of the QCD Lagrangian to include two spurion fields qv, QA H] 

£q™ = $lY% + QL F^ L + + q R F^ R - 4>(s - i^p)^ , (1) 

with 

lfi, = Vp- a M , q L = qv ~ qA, r ll = v ll + a M , q R = q v + qA ■ (2) 

The number of light flavours is assumed to be N F = 3, and the sources in eq.@ are 
traceless 3x3 matrices. Originally [I], the spurions qv, qA were introduced only for the 
purpose of classifying the independent terms in the chiral Lagrangian and afterwards one 
would set qA = 0, qy = Q. We will exploit here more general applications by consid- 
ering the spurions on the same footing as the ordinary external sources. We note that 
renormalization will generate terms which are non-linear in the sources. Those which are 
quadratic in the spurion sources are of particular interest for our purposes and will be 
discussed in detail in sec0J It will also be argued that no difficulty arises in the definition 
of the generating functional in the presence of the qy, qA spurions up to quadratic order 
in the spurions. We now introduce a set of five correlation functions defined by taking 
two functional derivatives with respect to u« or and two derivatives with respect to the 
charge spurions qy or qA- The first one is defined as follows 

where W is the generating functional of connected Green's functions. The remaining four 

2) < A a a A b (3 Q v Q d v > 



4 



3) < V^Q A Q d A > 

4) < V°V$Q c v Q d v > 

5) < A a a V h p Q A Q d y > (4) 

are defined in an exactly analogous way. These correlators involve two Lorentz structures 
g a p and k a kp. For our purposes, we can restrict ourselves to the Lorentz structure propor- 
tional to g a R and we will consider the limit k — > 0. Before taking this limit, however, we 
must consider two possible kinds of singularities: a) l/k 2 pion poles and b) log(fc 2 ) cuts 
due to pion loops. The latter singularities drop out in the leading large N c approximation, 
which we will frequently invoke. A pion pole can appear only in one flavour structure 
of the fifth correlator < A^V^Q A Qy >. Apart from this particular case, we can simply 
set k = 0, meaning that we restrict ourselves to constant vector and axial-vector sources. 
This leads to great simplification since all terms like duV^, d^a v etc... can be dropped in 
the chiral resonance Lagrangian. Similarly, the qy and q A sources can be taken as constant 
for the present applications. 

From the definition, eq.(J2J) and the explicit expression of the terms linear in the sources 
eq. Q we can easily see that the correlators introduced above can be represented as a convo- 
lution of an ordinary connected QCD 4-point Green's function and the photon propagator. 
For example, one finds 

< A a a A b pQ A Q d A > = J d 4 xd A yd A zexp(iky) < 0\T(A a a (x)A b p (y)A c p (z)A d a (0))\0 > D p /(z) . 

(5) 

Because of the short distance singularities, the integral over d A z is expected to be divergent. 
The divergence is removed upon including the additional contributions stemming from the 
counterterms and one may use d-dimensional integration as a regularization method. 

Let us now examine the flavour structure. The correlators depend on four flavour 
indices. At order e 2 p 2 one can show by using the chiral Lagrangian that there are only 
four independent tensor structures. It will prove convenient to use the following basis 

jtac jtbd _|_ jtad j^> c 

e 2 = 5 ac 5 bd + 5 ad 5 bc 
e 3 = d tac d tbd + d tad d thc 

e 4 = fabftcd _ ( 6 ) 

Other tensor structures that arise in practical calculations can be expressed in terms of 
this basis by utilizing SU(3) trace identities. The following relations are useful 

d tab d tcd = 1 1 _ 1 

6 3 2 A 
6 ab 6 cd = \b x + . (7) 

We can now expand the correlators over the flavour basis. The first four correlators are 
expanded over ei, e2, e3 because of the permutation symmetry a <-> b, c <-» d. For 
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instance, one has 

3 

< A a a A\Q c A Q A > = ig aP Y, < AAQ A Q A > jej . (8) 

3=1 

The expansion of the last correlator involves also, 

4 

< A%V b p Q A Q d v > = ig aP Y < AVQ A Q V > jej . (9) 

3=1 

This expansion apparently generates sixteen independent flavour components. As we will 
see below there are only ten components which appear at the chiral order e 2 p 2 considered 
here. 

2.2 Correlators at order e 2 p 2 

We start by computing the correlators defined above©® in the chiral expansion up to 
the order e 2 p 2 . Setting the sources s and p equal to zero, the leading order Lagrangian 
collects the terms of order e 2 and p 2 

\F, V F^ - ^(d,Fn 2 + \M 2 F^ 

(10) 
with 

D tl U = d„U-i(r fl + q R F l ,)U + iU(l tl + q L F^ . (11) 

The photon is endowed with a small mass (which may be counted as 0(p)) in order to 
regulate infrared divergences which can appear in the chiral limit. We will also need the 
following part of the chiral Lagrangian of order e 2 p 2 [3] 

= \KiF 2 (D^UD^) (q L q L + q R q R ) (12) 

+K 2 F 2 (D fl UD^) {q L tfq R U) 
+K 3 F 2 ((D f ,Uq L U^) 2 + (D^WqRU) 2 ) 
+K 4 F 2 (D^Uq L U^(D^q R U) 
+K 5 F 2 (D^D»Uq L q L + D^UD»Wq R q R ) 
+K 6 F 2 (D ll U^D»Uq L UURU + D D*tf q R U q L rf) 
+K 12 F 2 (UD^[V^q R ,q R ] + D^U[V ^q L ,q L }) (13) 
+K 13 F 2 (V ll q R UV fl q L 17+) 
+K li F§ (V^ L V^ L + V^qaWqii) 

where 

V M gL = d^qL - q L ], V ^q R = d^,q R - q R ] . (14) 



,(2) 

"chir 



C" 



i re 

chir ' chir 



Cfhir = C{q L tfq R U) , 



(3 



(2) 

As usual, the correlators collect the contributions at one-loop generated from £V- an d 
the tree contributions from jC e cf ^ r . Let us first quote the results, some details about how 
the calculation proceeds will be given later. For the components of < AAQaQa > one 
obtains 

< AAQ A Q A > 1= F 2 [2(K{ - K r 2 ) + 2(K r 5 - K r 6 ) + 4K[ 2 + K{ 3 + 2K r u + Xl (M 7 , M ) 

< AAQ A Q A > 2 = F 2 [-2(2^ + KI) + - + ^(M)] 

< AAQ A Q A > 3 = F 2 [6(K[ - Kl) + 2(F£ - XJ) + -Z% //)] (15) 

where KJ stand for the renormalized scale dependent parameters. The function X i(M 7 , /i) 
is one of Xii * = 1)2,3,4 which are the one- loop chiral contributions involving a virtual 
photon. We will describe the computation in detail in the next section. The contributions 
proportional to the function Z(k, fi) arise from purely pionic one-loop diagrams having 
one vertex from the term proportional to C in ea.(|lU|), This function reads 

Z (*'"'=16^(- l0 ^ + 1 ) ' < 16 > 

While we include this contribution here for completeness we note that it is of subleading 
order in the large N c expansion. In the following, we will drop such contributions for 
consistency with the leading large N c approximation which we will make in the treatment 
of the resonances. We also note that the correlators must be independent of the chiral 
renormalization scale /i. This constraint allows one to recover the scale dependence of 
KKfj,) computed in refs.|i] [37]. generalized to an arbitrary £ gauge. The components of 
the second correlator has the following expression 

< AAQ V Q V > 1= Fq [2(K{ + K\) + 2(K r 5 + K r 6 ) + AK{ 2 - K[ 3 + 2K r u + X2 (M T) M ) 

< AAQvQv >2= F 2 [~2(2K r 3 - K r A ) + ^(K r 5 + K r 6 ) - h(k, //)] 

< AAQ V Q V > 3 = Fq [6(K[ + K r 2 ) + 2(K r 5 + K r 6 ) - ^Z(k, fi)} . (17) 
For the next two correlators, only the first flavour structure appears at this chiral order 

< VVQaQa >x= F 2 \-K{ 3 + 2K r u + xs(M 7 , /i)] 

< VVQvQy > 1= Fq [K{ 3 + 2K r 14 ] . (18) 

The last correlator has two non-vanishing flavour structures. The structure proportional 

(2) 

to e4 is the only one to receive a tree- level contribution from C cl ^ ir , 

< AVQaQv >i= F 2 [2^[ 2 + + X4(M 7 , /x)] 

2(7 

< AVQaQv >4= -p- + F 2 [-F[ 3 + 2F[ 2 + X4 (M 7 , //)] . (19) 
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We now introduce a set of combinations of these correlator components which enjoy 
useful properties 



III = 


< AAQ A Q A 


>i - 


2 < AVQaQv 


>i + < VVQaQa >i 


n 2 = 


< AAQ V Q V 


>i- 


2 < AVQaQv 


>! + < VVQvQv >i 


n 3 = 


< AAQ A Q A 


>3 






n 4 = 


< AAQvQv 


>3 






n 5 = 


< AAQaQa 


>2 " 


2 

- < AAQaQa 


>3 


n 6 = 


< AAQvQv 


>2" 


2 

- < AAQvQv 


>3 • 



These six combinations are in one to one correspondence with the parameters K^,...,Kq 
which we intend to determine. One can show (see that they are independent of the 
gauge parameter £, and we will show below that they get no contributions from the QCD 
counterterms. This last fact implies that the photon loop integral must be converging. 
Convergence must be understood in the sense that integration is to be performed in d 
dimensions and the limit d — > 4 is well defined. 
In the leading large iV c approximation one has 

n 5 = n 6 = (large N c ) . (21) 

This is not difficult to show. The graphs which are of leading order in N c are the planar 
ones which involve a single fermion loopjHHj- This implies that they involve a single trace 
over the flavour A matrices. Furthermore, we are interested in four-point QCD correlators 
which are symmetric under permutation of the Lorentz indices a, (3 as well as in the other 
two indices pa (see eq.lfKjl). They must then have the following structure 

< A a a A b p Q A Q d A > = ig a(3 (E{{X a , A b }{A c , X d }) + F(X a X c X b X d + A a A d A fe A c >) 

= i 9aP (2(2E - F) ei + ^(2E + F)(e 2 + |e 3 )) (22) 

making use of relations (J7J), where E and F are arbitrary, and the same holds for the 
propagator with QaQa replaced by QvQv- This proves relations (J2Tj) . From these, one 
can recover the large N c relations among the corresponding K\ parameters 

K\ = -K{ K\ = 2K r 2 (large N c ) (23) 

noted in refs.jl] and [30] • We will further consider the following three combinations 

n 7 = <VVQ v Qv >i- <VVQ a Qa>i 
n 8 = < AVQaQv >i- < VVQvQv >i 
n 9 = lim k 2 < AVQaQv > 4 (24) 

From these combinations we will recover established results concerning the couplings C, 
1^12) K[ 3 . Hi and Ilg depend on the gauge. Convergence of the photon loop integral can 
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be established for II7 and Ilg while for Ilg it holds in the particular gauge £ = and in a 
limited sense, as we will discuss in more detail. This limitation affects the determination 
of K\ 2 - We will not attempt to estimate the coupling if[ 4 which is a pure source term. 

One further remark is in order. One usually expects the coupling constants in ChPT 
to be associated with QCD correlators which are order parameters, except for the pure 
source terms. Looking at eas. (|20j) this, at first sight, seems sometimes not to be the case 
here. For instance, II3 is associated with the correlator < A^A^A^A^ > which does not 
vanish in perturbation theory. However, this is simply due to the fact that the components 
< AVQaQv >3 and < VVQvQv >3 vanish at chiral order e 2 p 2 , so one could equally 
well express II3 in terms of a correlator which does behave as an order parameter. 



2.3 Calculation of the one-loop photon contributions 

In this section, we will compute contributions to the generating functional W(v^, a^, qy, qA) 
arising from the chiral Lagrangian C c ^ ir . The use of the spurions qy, qA as sources gives 
rise to a simple diagrammatic method which we will in this section and the following 
ones. As there were some errors reported in related calculations in the past (see e.g.[36j) 
and we will, in fact, find some discrepancies with BU, we will present our calculations 
in some detail. Here, we will compute the functions Xi{M~,fi) which appear in eqs. (|17j) . 
They arise from sets of one-loop diagrams containing one photon line and one pion line, 
the vertices being obtained from the leading order Lagrangian IjlOj) . The corresponding 
standard £-gauge massive photon propagator reads 

^w-5^){^ + «-i)p^j}- < 25 > 

The list of the vertices which are effectively needed is displayed in the appendix. We 
remark here that these contributions are to be included both in the chiral expansion of 
the correlator as well as in the expression from the resonance Lagrangian since the latter 
includes the piece from eci. dTUl) . 

Let us now compute the contribution to < AAQaQa >%■ A single diagram contributes, 
shown below 



The following conventions are used to represent the sources 
• = qA source -)- 



o = qy source x = v ii source 

The relevant vertices are listed in the appendix. One must keep in mind, for instance, that 
the coupling of one axial source to the pion does not contribute to our correlators. The 
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photon line which joins the two q-sources is not explicitly drawn. Each diagram is easy 
to interpret in terms of the integral representation eq. (JSJ) . We will first give the results in 
integral form, in order to display the integrands, in which we perform the replacement of 
PaPfS by p 2 g a /3/d. In integral form, this diagrams gives 



< AAQaQa >i 



—i dp 
(2vr) d p 2 (p 2 -£M 2 ) 



(26) 



and the other components < AAQaQa >j get no contributions. Next, we consider 
< VVQvQv >i and, again, there is a single diagram to compute shown below 



which gives 

< AAQ V Q V > 1 



1 



-id d p J 1 

~{2^ 1 p 2 (p 2 -M 2 ) ' d (p 2 - M 2 )(p 2 - CM 2 ) 



(27) 



and the other components vanish. One also finds that < VVQyQy >i receive no contri- 
bution from pions. Next we consider < VVQaQa >i, in this case there are four diagrams 



X 



X 



X 



Calculating these diagrams, one finds that the last two cancel each other and the first two 
give the following integral 



< VVQaQa >i= F 1 



— 

"(2<F" (p 2 -M 2 )(p 2 -eM 2 ) ' 



(28) 



'7/^ s 7' 

Finally, we have to compute < AVQaQv >i- There are four diagrams in this case 
« X» xf- X-f- 



c 



-Ho 



•o 
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Computing these diagrams gives 



< AVQ A Q V > 1= -F 2 



-i dp J 1 



1 



1 



(2vr) d 
1 



d p 2 (p 2 - £M 2 ) 2 p 2 (p 2 - M2) 



< ^FQaQv >4= 



2d{p 2 -M 2 ){p 2 

F 2 (l - d) y -id d p 1 
jfe2 y (2vr) d f- 



+ < AFQaQv >i 



(29) 



The 1/A; 2 term is generated by the third and fourth diagrams and the integral vanishes 
in dimensional regularization. We can now compute the d-dimensional integrals and the 
divergences are removed using the chiral MS prescription [3] • The results for the loop 
functions Xj(M 7 ,/x) are collected below 



Xi (My,//) 

X4(Af 7 ,/i) 



1 



167T 2 

1 

167T 2 V 

1 



M 2 

£log^ + £log£ 



,2 



2(£ + 3)log^ + 2£log£-£ + l 



16vr 2 
1 

167T 2 



,2 



6(£-l)log^ + 6£log£ + £-l 



M 2 

2(£-3)log^ + 2£log£-£-l 



(30) 



Using these results and those from the preceding section which display the contributions 
involving the chiral couplings K\ we can deduce the scale and gauge dependence of the 
latter. One finds that K{,...,Kq are independent of the gauge parameter £ and the scale 
dependence is in agreement with the results of 01 EZ1 - The couplings K[ 2 , K[ 3 , K[ 4 do 
depend on £. Defining the beta functions Sj in the conventional manner 



dfi 



1 



16vr 2 



one finds from eas.(|3U|) 



:i2 = j(3-2f), E 1S 



-11 



•s(i-0- 



(31) 



(32) 



These results agree with those of Urech in the special case £ = 1, and they reproduce the 
results presented in ref.[SJ where a different set of correlation functions were employed. 



3 Correlators from the minimal resonance Lagrangian 

3.1 0(p 4 ) resonance Lagrangian 

Our starting point for determining the K\ parameters is the integral representation eq.© 
in terms of 4-point correlation functions. The first approximation that we make is to 
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consider the leading large N c limit. In this limit, it is an exact statement that QCD 
correlators can be expressed as a sum over tree graphs involving resonance propagators [SSI 
(see e.g.jSH] for a review of recent applications). The second approximation that we make is 
to restrict ourselves to contributions from a finite number of resonances. This is acceptable 
in our case because we consider combinations of correlators which fall off as powers for 
large momenta. 

In order to guarantee the validity of chiral Ward identities one may compute the 
correlators starting from a chiral Lagrangian with resonances. In cases where the integrals 
can be shown to converge, imposing the convergence property generates a set of constraints 
for the parameters of the resonance Lagrangian which we will call consistency conditions. 
These conditions will be discussed in detail in the following. We begin with the resonance 
Lagrangian which was used by Baur and Urech[2l?] 

£res = £-chir + ^S,S + ^V,A ■ $3) 

BU have used exactly the Lagrangian proposed by Ecker et al. [34] (except for leaving out 
the 7r' resonance multiplet) which we will call the minimal resonance Lagrangian. Let 
us recall the notation and formalism used in ref. |34j . The formalism that is used is one 
in which all resonances transform homogeneously under the non linear representation of 
the chiral group including the vector and the axial-vector resonances which are described 
as antisymmetric tensor fields. In this approach, all the resonance fields can be assigned 
a chiral order equal to two, such that the resonance Lagrangian of ref. [31], contains all 
possible terms, for each resonance field, which are of chiral order four. The only terms of 
chiral order six which are present are the kinetic energy terms. 

We recall here the details: in the scalar sector an octet S and a singlet So are considered 
and the Lagrangian involves two coupling-constants, Cd and Cd, 



±(V x SV x S-M 2 S 2 ) + ±\ 
+c d (Su^) + c d S (u^) . (34) 



£s,s = ^(V A 5V A 5 - Mis 2 ) + -[d x s d x s - M 2 So S 2 ] 



Next, a nonet of vector and axial- vector resonances are considered and the minimal La- 
grangian involves three coupling constants Fa, Fy and Gy, 

£va = -\ E (v A i? AM v^- \m 2 r r^r^) 

1 R=V,A 1 
F A , A ?ix V \ , Fy : iGy 



'2 (V/ - } + 27i (wr) + 27i (VKX]) • (35) 

We will compute the basic correlators introduced in sec.l first from this minimal resonance 
Lagrangian. It is convenient to keep using the qy, qA source spurions at this level, intro- 
duced in association with the and sources, with the replacements 1^ — > + qiF^, 
r n —> ?> + QrF[i- The diagrammatic calculation is exactly equivalent to using the convo- 
lution representation eq.© in which the QCD correlator is computed from the resonance 
Lagrangian. In principle, we should be able to reproduce the results of BU (who used 
a completely different approach to determine the K?'s ). We follow their convention to 
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regularize the UV divergences in the chiral MS prescription and we introduce a regu- 
larization scale A. The various vertices needed in the calculation are collected in the 
appendix. We will use the photon propagator in a general £ gauge as in eq.(|25|). Since 
the resonances are massive one can set the photon mass M 7 = in the following. Since 
the 0(p 2 ) chiral Lagrangian is also part of the resonance Lagrangian, the pion contribu- 
tions Xi(My, A), i = 1,2,3,4 which we already computed (eas.(|3Uj)) must be added to the 
resonance contributions. 



3.2 Contributions proportional to F\ and F v 

We will now consider successively the contributions proportional to F\, Fy, GyFy, Gy, 
GyFy which are the only ones generated from the vector and axial- vector part of the 
resonance Lagrangian ea. (|33|) . One easily sees, at first, that the following correlator com- 
ponents < AAQaQa >i and < VVQvQv >i get no contributions proportional to F\. 
Next, one finds that there are five diagrams which contribute to < VVQaQa >i, the 
relevant vertices can be found in appendix A. The diagram shown below is found to make 
a vanishing contribution 

F A • 



F A 



as it is antisymmetric in the a (5 indices. The remaining four diagrams are 
F A • 



The result is of the following form 



< VVQaQa >i = F aJ E A(P) , (36) 



with the individual diagrams contributions to the integrand reading 

1 s / 2-d 6-d 4M| 



D l = i 1 ~ 3) T7273 ^72T + 7Z3. + 



d'\Ml{p 2 -M\) {p 2 -M\Y (p 2 -M]) 3 



1 -2 
2 " (1 " d\p*-M\y 
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2-d 



+ 



d> \M 2 A {p 2 - MV) p 2 {p 2 - M\) 



Da = (1- 3) 



1 / 2(d-2) 



d ' \ M\{p 2 - Ml) (p 2 - M\) 2 



(37) 



respectively. Next, we consider < AAQyQy >. A single diagram contributes (with a 
crossed one) which gives the same result as D3 above. Finally, we turn to < AVQaQv > • 
There are three diagrams which contribute (and no crossed diagrams in this case), which 
are shown below. 



F A • 



+ 



F A 



+0 



: 



These diagrams give 

< AVQaQv >i = F\(l 
k 2 < AVQaQv > 4 = F\ 



—i dp 
(2ir) d 



F A d 

F A 

d + Z-1 1 

2p 2 {p 2 -M 2 A ) {p 2 -M\) 2 



-i d p 1 — d 



(27r) d p 2 -M\ ' 

Computing all the integrals, in summary, the following results are obtained 



< AAQaQa >T a = 



< AAQyQy > x 



3 F 2 



F 2 



< VVQaQa >i 

< VVQvQy >[ 

< AVQaQv >T a 



416tt 2 
3 Fl I 



F 2 
1 A 



2167T 2 





M jt 
A 2 " 



3 F 



1 



M 2 A 



8167T 2 

M 2 F 2 2 



k* < AVQaQv >? = 3^(| + log ^) + 0(A; 2 ) . 

lo7r z 3 A z 



(38) 



(39) 



There is a simple correspondence between the diagrams proportional to Fy and those pro- 
portional to F\ considered above so we give directly the results for the various components 
of the correlators 



< AAQaQa >i 



F 2 



3 Fl 



4167T 2 



(1 + (i + b, 



A 2 
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< AAQ V Q V >, 



< VVQ A Q A >; 



1 v 



< VVQvQv >\ 



< AVQaQv >i 



k 2 < AVQaQv > x 






3 F 2 



2 16vr 2 



3 Fj 



v 



8 16vr 2 



(1 + U+lo 



M 2 \ 



A 2 



F' 2 
i 



M£i<£ 2 My- , 



3.3 GyFy contributions 

We have two diagrams which contribute to < AAQaQa >i shown below. 

Gy 

G v 



(40) 



This results in the following integral for the correlator 



oWQaQa >i = 4GyFy(l 



1, 



i d d p 



d-2 



(2vr) c 



M 2 (p 2 



m 2 ; 



+ 



1-d 



p 2 (p 2 - M{ 



(41) 



Next we have one diagram contributing to < AAQvQv >i which is non vanishing 

-ho 




and which gives 

< AAQvQv >, = ACyFy I =0 vH /_ Mi) ■ (42) 

The two correlators < VVQaQa >i and < VVQvQv >i are then found to get no 
contributions proportional to GyFy. 

Next, we consider < AVQaQv >%■ One first finds a set of 7 diagrams which contribute 
both to < AVQaQv >i and < AVQaQv >a- The three diagrams drawn below give zero 
contribution 
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G v 



F v 

The next four diagrams are 
Gv^P 



4# 




Computing these four diagrams yields the following result for < AVQaQv >i 



< AVQaQv >i 



G v F v (l-~: 



—id d p I 



-d 



My] 



+ 



(P 2 



Myf 



(43) 



these diagrams make an identical contribution also to < AVQaQv >&■ I n order to get 
the full contribution to < AVQaQv >a we must also add the following two diagrams 

HlX -IK 



Gv 



Fv 

These diagrams are slightly more complicated to compute because one must allow for a 
small momentum to flow through the vp source and expand the result in powers of k 2 . 
One finds that the k 2 pole vanishes and then adding the constant contribution to that 
arising from the preceding four diagrams one finds that they cancel. In summary, the 
contributions proportional to GvFv which are non vanishing read 



< AAQaQa >i 



< AAQvQv >i 



GyF\/ 



SGvFvfl Mf\ 

16n 2 \6 8 A 2 J 

_3_GvFv (I Mf 

16^ 2 l 6 g A 2 



2 16vr 2 \6 



A 2 



(44) 
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3.4 



Gy and G V F V contributions 



These contributions are simple because a single correlator component is found to be non 
vanishing: < AAQaQa >i in the case of G v and < AAQyQv >i hi the case of G V F V . 
There are three diagrams proportional to Gy 



G v 



Gy 



G v 



G v 



G v 
G v 



: 



and a single one proportional to G V F V 



F v 

G v 
G v 

F v 



After a small computation the following results are obtained, 



< AAQ A Q A >i= 4G 2 V (1 



', d d p 



d 



{2ir) d \M v (p 2 



My] 



< AAQyQv >r- 



{ G V F V 
Fo 2 



(1-3) 



+ 



d- 1 



p 2 (p 2 



My) 



—i d d p 1 
~(2^y r {p 2 - M 2 ) 2 " 



In summary, the non vanishing contributions proportional to Gy read 



(45) 



< AAQ A Q A > 1 



Gy 



3Gy(l^ ] My' 

mT 2 U + log 



< AAQyQv >x 



GyF v 



3G 2 V F 2 1, M 2 



16vr 2 F 2 \ 6 



n + l0 § 



A 2 



(46) 



3.5 Contributions proportional to cj, Cd 

The scalar singlet turns out to contribute only to < AAQaQa >2 an d the scalar octet to 
< AAQaQa >3- The three Feynman diagrams to be computed are depicted below. 



1 
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The dependence on the photon gauge parameter £ is found to cancel out, and the following 
results are obtained for the singlet and octet scalar contributions 



< AAQaQa > 2 = 1«3 (l - I) / =0 ^mI-^ 

< AAQ A Q A > 3 = 8cl (l - i) / ^ ^^-^ . (47) 

As before, we have replaced the term p a P/3 in the integrands by p 2 g a p/d. There are three 
additional tadpole diagrams. Let us quote the result for completeness 

, A a A b n c n d ^tad_ ■ ( & c d jsab jscd , ^^d x ab x cd\ f ~jj P ^ f A a\ 

<A a A p Q A Q A > - l9ap y—d d + ] ^5 5jy_ 3 -_. (48) 

Computing the integrals we obtain the following results for the scalar resonance contribu- 
tions which are non vanishing 

c 2 , ( M% l\ 
<AAQ A Q A >f al = l2j^ [log^ + gj 

r 2 ( /If 2 1 \ 
< AAQ A Q A >r l = 6^ (fog + (49) 

One notices that no dependence upon the gauge parameter £ appears except in some of 
the terms proportional to F A and -Fy. This is to be expected because, quite generally, the 
part proportional to £ in the integral representation can be simplified by making use of 
chiral Ward identities and expressed in terms of 2-point Green's functions |31j. 



3.6 Comparison with the results of Baur and Urech 

Up to this point we have used exactly the same resonance Lagrangian terms as Baur and 
Urech|2*9"] and we are therefore in a position to check their results for the chiral parameters 
K\. Making use of the relations (|15j) - (|19ft . for each resonance term contribution to the 
correlators (j3J)@ we can solve for the -KJ's which will then depend on the two scales ji 
and A. In ref.j22] the values of the two scales have been taken to be equal. Following this 
convention for the sake of comparison, and using the results derived above, we reproduce 
the results of ref.[2H] for the parameters -ff -Kg. 

We disagree, however, on a number of contributions which concern KI2, K [ 3 , K\^. 
Contrary to BU, we find that the following quantities are vanishing. 

„ „ „ „ (-il /~t2 ri2 (~i2 t?2 (~i2 t?2 ri2 t?2 

-"-13 — ""-14 — -"-12 — -"-13 — -"14 — -"12 — -"13 _ -"l4 ~ u ■ V 0U J 

In the case of if [3, K[ 4 this follows from where it is shown that they can be obtained 
from 2-point correlators < QvQv > an d < QaQa >• Obviously, these can only pick up 
contributions proportional to F A and F v , which is indeed what we recover here. Similarly, 
one can also obtain K[ 2 from a correlator with three currents < PQaQv > which picks 
up no contribution proportional to Gy. 
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4 Short distance constraints 



We now need to determine the divergence structure in the integral representation of eq. JSJ 
according to QCD. One way to proceed is to determine the necessary counterterms which 
are generated in QCD with a source term as in eq.([T)). where is a dynamical (massive) 
photon. This QED-like sector differs from ordinary QED by the fact that we have replaced 
the charge matrix Q by a vector (qv) as well as an axial-vector {q^) spurion source. We 
must consider this theory as a kind of effective theory and we assume that the divergences 
can be absorbed into a set of local Lagrangian terms. We will restrict ourselves to an 
expansion quadratic in the spurion fields. This represents an important simplification 
since, for instance, we will only need the free photon propagator and thereby avoid the 
problems of the photon propagator at higher order in the presence of axial couplings |39j. 
Another restriction is that we will only consider contributions which are of leading order 
in N c . Dimensional regularization together with the naive 75 prescription will be used. 

4.1 QCD+QED counterterms 



(b) 




Figure 1: Typical divergent diagrams quadratic in the spurion sources (a) with an open 
quark line (b) with a closed quark line. 

Let us consider QCD in the presence of dynamical photons and a source term of the 
form of eq.(^Q). The effective action is expanded to quadratic order in the spurions qy 
and qA and we drop the contributions which are subleading in N c . We can then split 
the divergent terms in the effective action into two categories, (a) one which collects the 
graphs with an open quark line and, (b) one which collects the graphs with a closed quark 
loop (see fig. 1). The former ones are those which will prove to be of relevance to us. 

Let us first consider the contributions in category (b). The graphs in this category, at 
any order in a s (in QCD the three gluon vertex has dimension p, the four gluon vertex has 
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dimension zero), have dimensionality p 4 . The corresponding local counterterms therefore 
must involve four derivatives. Examples of such terms are, 

(VVVV + VVVV), G^GHVpfcVV + V^V^), etc... (51) 

If we restrict ourselves to the correlator combinations H , eq. Q2UJI (from which the contri- 
bution proportional to K14 drops out) then all the contributions from such counterterms 
also drop out or are suppressed in the large N c counting. 

We are finally left with the contributions from the first category. Taking into account 
the structure of QCD, all graphs in this category have dimension p, implying that the 
counterterms must involve one derivative. Furthermore, they must involve two spurions 
of the same kind (either qi or q R ). One can form only two such independent terms 

% — — 
£ = ^ z 2 {4>L[qL, V M g i ]7 At ^ i + TpR[qR,V fl qRh^R} 

+\Z' 2 {^ql^D^L + i> R q 2 Rltl D^ R + h.c) (52) 

where 

iD^ L = {id" + 1^ + G^l, iD^r = {id" + r„ + G^ R . (53) 
At order zero in the strong coupling constant one has the well known QED result 

Z 2 = Z' 2 = ^(-e)[r(e) + log(47r)] + 0(a s ) . (54) 

Equality between Z 2 and Z 2 need not hold beyond the lowest order. 

Let us now discuss the contributions proportional to Z 2 and Z 2 to our basic correlators 
(jHJ) © • The first result is that the contribution proportional to Z 2 vanishes at leading order 
in N c . This can be shown in a pedestrian way. Because of the conservation of chirality 
in one fermion loop graphs, the contributions concern only the correlators < LLQlQl > 
and < RRQrQr >, it suffices to consider the former. For a given gluon field configuration 
G fJi (x) we can express the functional derivatives with respect to the sources in terms of 
exact fermion propagators Sg{x,u) and we obtain the following expression 

< L^QlQt >G = ^J d *xd i y' 

-S(y) < X b X a {X c ,X d } > < 7 a PLS G (x,0)^PLS G (0,x) > 

-6(x) < X a X b {X c , X d } > < 7 ' 3 P L S G ( ? /,0)7 Q P L 5 G (0,2/) > 

+ < X a X b {X c ,X d } > < 7 Q P L S G (x,y) 7 /3 PL5 G (y,0)[^+^]P L 5 G (0,x) > 

+ < X b X a {X c ,X d } > < 7^P L 5 G (y,x) 7 Q P L 5 G (x,O)[^+0']P L 5 G (O,y) > ] , (55) 

where Pl is the projector Pl = 1/2(1 — 75) and the notation p means that the derivatives 
are acting on the right. There is another similar contribution coming from the hermitian 
conjugate part of the Lagrangian. Using the equation obeyed by the Dirac propagator 

[0+ G]P L S G (x, y) = P R 5\x - y) , (56) 
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we indeed find that the various terms in ea. (|55|) cancel each other. We have thus shown 
that the only relevant contributions from the counterterms are those proportional to Z 2 . 
The calculation of these is performed by computing the functional derivatives according 
to the definition of eq.(j3J) and using the following result for the two-point correlators that 
appear 

J d 4 xe ikx < A s Q (x)4(0) >= iF$6 st - g a p\ , J d 4 x < V°(x)V^0) >= . (57) 

We find that the counterterms contribute to the following components of our basic corre- 
lators 

< AAQ A Q A >?= -Z 2 Fl 

< AAQ V Q V >f = -Z 2 Fl 

< AVQ A Q V >?= -\ziFl 

< AVQ A Q V >t= - X -ZiFl . (58) 
There are no contributions to the other components. 



4.2 Summary of the log A 2 dependence 

Let us consider the combinations of correlators II j introduced in eas.(|2U |) (|24 |) . In all these 
combinations the contribution from the coupling drops out. In the combinations 
IIi,...,n6 as well as II7, II9, the contributions from the renormalization parameter Z 2 also 
drops out, as is seen from using eq.fJSSJ). As a consequence, these quantities must be finite. 
More precisely, finiteness holds in the sense of dimensional regularization which was used 
to establish the form of the counterterms in the preceding section: one must compute the 
photon loop integral in d dimensions and then take the limit d — > 4 which must be finite. 
From these eight combinations one can then determine the chiral couplings K[,...,Kq as 
well as K[ s and C. The situation is different for lis since, in that contribution 
proportional to Z 2 is present. The integral is not finite in general. It is finite to zeroth 
order in the strong coupling a s when the gauge parameter £ = 0. This limitation affects 
the determination of the coupling K[ 2 ■ 

Let us examine how these finiteness conditions are satisfied when the relevant corre- 
lators are computed from the minimal resonance Lagrangian of eqs. (|M4() . (jH5jl . Collecting 
the results presented in detail in secOH we obtain the following dependence upon the 
regularization scale A 



Tl l ~— 2 -^- T [FZ + Fl-2Gl] 

U2 ~- 3 2 [F ° 2 + ^ - 4GvFv + 2G ^f ] 

2 



rr « lo £ A r 21 



n 



4 
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TT ^Og^ro^ 2i 

5 ~ 16^ [3 ° d ~ Cd] 

n 6 = o 

3 log A 2 2 2 2 

7 4l6^" [ v ~ aJ 



3logA 2 r ^ 2 p2 P 2 

J ^^r^ M v F v-MlF\] . 



At this level, there is at first one consistency condition which is satisfied, but trivially so, 
since II4 is identically vanishing. In the sequel (sec !5.6|) we will consider physically relevant 
interactions and attempt to satisfy the condition in a less trivial way. The consistency 
conditions for II7 and Ilg can be satisfied and they generate the two Weinberg sum rule 
relations 

Fy = Fq + Fl M V F V = M\F\ . (60) 

These relations will be assumed to hold from now on. Furthermore, we will assume the 
following relation to hold, 

GyFy = Fq (61) 

which follows [53] from imposing that the asymptotic behaviour of the vector form- factor 
of the pion 

lim tF%(t) = (62) 

t— >oo 

is satisfied. Using this, we observe that there are only two consistency conditions which 
are fulfilled, which are those associated with II7 and H9. The corresponding evaluation of 
II7 and Ilg gives 



n 7 



3 F 



4167T 2 



, M v ^ 

log i^ + 2 -^T 



where z denotes the mass ratio 

* = ^ (64) 

My ^ ' 

which one expects to be approximately equal to 2. We expect 1I5 and H6 to vanish as a 
result of the leading large N c approximation. This condition is trivially satisfied for Tl§ 
and it can be imposed for II5 by setting 3c 2 , = c 2 . 

5 Beyond the 0(p A ) resonance Lagrangian 
5.1 Survey of 0(p 6 ) terms 

In order to comply with the other consistency conditions, it is necessary to enlarge the 
resonance Lagrangian of ea. ()33|) . Since this Lagrangian (for a given resonance content) 
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contains all terms of order p 4 , it is natural to investigate terms of order six. Generically, 
let us designate a resonance building block by R (which counts as p 2 ), a chiral building 
block by tt (which counts as p) and a source term by J (which counts as p 2 ). In this 
schematic notation, there are four classes of terms of order p 4 which involve at least one 
resonance 

RR, RJ, Rtt, Rtttt . (65) 

At order one can have the same terms but involving two more derivatives and, in 
addition, one can have terms formed from three, four as well as five building blocks 

RRR, RRir, RRJ, RJJ, RJir, RRirir, RJtttt, Rtttttt, Rtttttttt . (66) 

We will not attempt to classify and discuss all possible remaining terms in this work. We 
will simply consider (except in sec !5.6l ) sets of terms of the form RRtt. The first set will 
be labeled as aipir and has the following form 

C aiP n = i(G ig ^g^ 3 + Chg^V 1 + G^ V 2 ) < [V w ^, \ 2 X K, > 
+i(H ig ^g^ + H 2 g^g^ + H 3 g^g^) < [V^, A x ]V^u^ > . (67) 



We will first show that by adding these terms to the minimal resonance Lagrangian we can 
reproduce the < VAP > 3-point function constructed in ref.[35 2 ■ I n addition, we will 
then show that it is possible to satisfy the consistency conditions associated with IIi and 
II2, and these provides a number of Weinberg-type sum rule relations. In an analogous way, 
we introduce an ao/ivr coupling (we assume nonet symmetry which reduces the number 
of couplings), 

C aofl7r = c A (S{V^A^,u u }) (68) 

which will allow us to satisfy the consistency condition on 1T3. 

Finally, we will consider a set of terms of the type Lvpn which contribute to H4. A 
calculation which is equivalent to an evaluation of this contribution was performed in 
ref.[Zj using the vector, instead of the anti-symmetric tensor, formulation for the vector 
fields. A classification of the terms in this latter formulation was recently presented by 
Ruiz-Femenia et al.|41j. we list below those of relevance to us, 



Cup* = W {ds({V x V^, ^ A K> + d 4 ({V°V^, V? x } 



C7 
My 



({V°V^,fi X }u x )} . (69) 



In enabled them to reproduce the 3-point function < VVP > as constructed in ref . |42j . 
The consistency condition on II4 is somewhat more complicated to satisfy in the anti- 
symmetric tensor than in the vector formalism. This will be explained in more detail 
below. 



2 An extension was recently considered |4()j which includes a tt' resonance multiplet and a corresponding 
resonance chiral Lagrangian was constructed. The result for the coupling constant combination G1 + G2 + 
2G3 which is needed in the following turns out to be unaffected by this extension. 
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5.2 C aipn constraints from < VAP > 



Let us examine the constraints on the six coupling constants of the aipn Lagrangian (|67|) 
which arise from the 3-point correlator < VAP >. This correlator is defined as follows 



W a p(p, q) 



d xd ye 



ipx+iqy 



<0\T(V^x)A 2 Jy)P%0))\0> . 



(70) 



Assuming an ansatz in terms of rational functions with a minimal number of resonance 
poles (i.e. one pion, one vector and one axial-vector resonance poles) it can be shown |3*T] 
that this correlator is uniquely determined from a) the chiral Ward identities, b) the 
constraint to match the leading asymptotic scaling behaviour, 



lim W a/3 {\p,\q) 

A— >oo 



with I 2 = (p + q) 2 



<uu> \{p a + 2q a )qp - q 2 g a p (p 2 - q 2 - l 2 )P a 



X 2 



n 2 



+ 



2p 2 q 2 l 2 



Qa/3 

p 2 q 2 l 2 



P a p = ppq a - p-qg a p 

Qap = p 2 q a q/3 + q 2 p a p/3 - p-qpaqp - p 2 q 2 9a/3, 



(71) 



(72) 



and, finally, c) the asymptotic constraint on the vector form factor of the pion ea. ()62j) . 
More detailed analysis of asymptotic constraints related to this 3-point correlator were 
performed in ref. |43| . The result, for < VAP >, reads 

(p a + 2q a )qp - q 2 g a p 



W a p(p,q) =<uu> 



q 2 l 2 



+ 



F{p\q\l 2 )P 



mi 



+ 



G(p 2 ,q 2 ,l 2 )Q a 



(p 2 -M 2 )(q 2 -M\)l 2 {p 2 - M 2 ){q 2 - M\)q 2 l 2 



(73) 



with 



F(p 2 ,q 2 ,l 2 ) = ^(p Z 



q 2 - I 2 ) + Mi 



V 



M\ 



G(p 2 ,q 2 ,l 2 ) 



-q 2 + 2M\ . 



Let us compute the VAP Green's function from the resonance Lagrangian C^ s +J0. ai fm- 
We find, for the F and G polynomials, 



(74) 



F(p 2 ,q 2 ,l 2 )=p 



+ q z 



n 



F, 



o 



F 2 - 2G V F V F V F A 



Fo 2 



G 2 - 2G 3 + H! + H 2 ) 



(G 2 + 2G 3 + Hi + H 2 ) 



1 

+ "^2 
^0 



-FyF A \ 



G(p 2 ,q 2 ,l 2 " 



[-G 2 - 2G 3 + H 1 + H 2 + AH 3 ] 

M 2 F\ - M\(F 2 - 2G V F V ) 
2G V F V 2F V F A 



F 2 



F 2 o 



-{-G 1 +H 1 + H 2 ) 



+ 2Mi^. (75) 
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Let us first compare the constant terms in eas. (|75|) and (|74|). Provided that the two 
Weinberg relations ea. (|6U|) are satisfied together with the relation GyFy = Fq, one finds 
that the constant terms are in exact agreement. Next, the following three independent 
relations between the coupling constants Gj, Hi 



Gi + G 2 + 2G 



H 1 + H 2 



Fa 
Fy 
F A Fy 



2F V 



2F A 



Gi 



2H, 



Fy 
2F A 



Gi 



(76) 



are found to be the necessary and sufficient conditions for the polynomials in eas.(|75[) and 
(|74j) to be completely identical. 

5.3 Linear contributions from C aifJ7T 

Let us now return to the 4-point correlators and consider the contributions which are 
linear in the couplings of ea. ()67|) . The couplings Hi do not contribute because derivatives 
of the sources are involved. There are contributions to < AAQaQa >i, < AAQyQy > x 
and < AVQaQv >i which we consider in turn. There is one kind of diagram relevant 
for < AAQaQa >i 

Fy, Gy+d 



G, 



J 



The upper vertex involves either Fy or Gy. These diagrams give 
< AAQaQa >i = -4F A (F V - 2G V )(1 - J =0^ 

+ 



{d-2)G 2 
M y {M\-p>) 



G 



123 



(M|- P 2)( M 2 -p2) 



where 



G 



123 



Gi + G 2 + 2G3 



(77) 



(78) 



denotes the combination which was determined from < V AP >. In the case of the 
correlator < AAQyQy > 1; there is again one kind of diagram which makes a linear 
contribution 
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G, 



J 



and the result is very similar to the preceding one 



< AAQ V Q V >! 



1 



-4^(1-3) 



: d d p 



d' J (2vr) 



+ 



G 



123 



(d - 2)G 1 

M\{Ml-y>) ' (Ml-p 2 )(M2-^ 



(79) 

Next, in order to obtain < AVQaQv > one must calculate a set of five diagrams which 
are depicted below, 

F A 



Ga 



4 



X 



Computing these diagrams gives 

< AVQaQv >i = -FaF v G 123 (1 



d' 



—id d p 
(2vr) d 



d-2 



+ 



2M\ 



+ 



(M 2 -p 2 )(M 2 -p2) 



2M£ 



(M^-p 2 ) 2 (M2 -p2) (M2-p 2 )(M^ -p 2 ) 2 
which is seen to be proportional to the coupling constant combination G123. 



(80) 



5.4 Quadratic contributions from C ai/ m and consistency conditions for 

iii, n 2 , n 8 

We turn now to the contributions which are quadratic in the couplings Gj. These concern 
only the two correlators < AAQaQa >x, < AAQyQy >i and they are given respectively 
from the two kinds of diagrams depicted below 
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The computation gives the following result for < AAQaQa >i, 



< AAQ A Q A >! = 4F|(1 - -) 



-i d p 



(Gi 23 )V 



+ 



(2ir) d \{Ml-p 2 ) 2 (M^-p 2 ) M v (M\-p 2 ) 



G\ (2-d)p* 

2\2 



(81) 

The result for < AAQyQy > 1 is obtained simply by replacing Fa by Fy and by inter- 
changing Ma, My and Gi, G2 in the expression above. 

We can now collect all the contributions which involve the coupling constants Gi, 



< AAQaQa >i 



3F A (F y - 2G V ) [ 

167T 2 \ 



-2G 2 z 



6 +1 ° g 



+ G 



123 



1 z log z , M T / 

— I 1- loe — — 

6 * - 1 8 A 2 



+ 



3^_ 
16vr 2 



4G| z - + loj 



Ml 
A 2 



(G 



123 J 



7z - 1 



log z 



6(z-l) (z-1) 2 



+ log 



A 2 



(82) 



Gi 



3F A F V f -2Gi /7 M 2 ' 



A 2 



+ G 



123 



Z lo£ 



(G 



123 J 



Z - 1 

7-z 



+ log 



M 2 , 
A 2 



+ 



3^ 
16vr 2 



,G 2 /5 



4-1 M v 
3 "A 2 " 



z 2 log z , M£ 

+ 7 iTTT+log^f 



6(1 -z) (z-1) 2 



< AVQaQv >i 



3 i ? A^VGi 2 3 
2 167T 2 



1 z log z 
- H — 

6 z- 1 



log 



K 

A 2 



(83) 
(84) 



Adding these contributions to those computed in sec|3] it is now possible to satisfy the 
consistency conditions associated with Hi, II2 and lis. Firstly one observes that the A 
scale dependence drops out from the combination Ilg automatically provided the first 
relation (|76j) holds for G123. Requiring that the consistency conditions for IT and H2 be 
satisfied yields two Weinberg-type relations for G\ and G2 

8G 2 - 4=Gi -1 = 



IG\ - -i=(2 - z)G 2 -1 = 0. 



(85) 



The previous determinations of Fy, Fa, Gy and of the combination G\ + G2 + 2G3 
have been used. Since these equations are quadratic, one could be concerned about the 
existence of real solutions. Real solutions do exist for any value of the mass ratio z and, 
a priori, one has a set of four different solutions. By studying the behaviour as a function 
of z we can reduce the multiplicity. We note that the Weinberg equations (|60|) have real 
solutions provided z lies in the range 1 < z < 00. Firstly, let z — ► 1: this corresponds to 
a situation where chiral symmetry gets restored and one must have G\ = G2 in this case. 
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This argument reduces the multiplicity to just two solutions, which we can write, 



4^ 
1 



(2- z + aJ(2- z) 2 + 2z), a = ±1 



(86) 



We can further reduce the multiplicity by observing the behaviour as z — > oo. In this limit, 
G\ remains finite while in the case of G2 only one solution remains finite corresponding 
to a = 1. Correspondingly, the divergence in III is logarithmic for a = 1 while it is linear 
for a = —1. The solution with a = 1, therefore, appears to be the most plausible choice. 
Using this solution, the following determinations are obtained for ITi, II2 and ITs 
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log 
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2z 2 + 7z 



10 



1)3 



2z + A 



n. 
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16vr 2 
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5.5 Contributions from £ ao / l7r and the consistency condition for II3 

By analogy with the preceding we now consider couplings which involve the scalar mesons 
of the type ao/ivr which will contribute to < AAQaQa >2 3- There is only one Lagrangian 
term which matters, shown in ea. ()68|) (assuming nonet symmetry). The term in which 
the derivative acts on u v does not contribute to our correlators nor to the physical decay 
amplitude which is discussed below. There are two kinds of diagrams which contribute to 
the relevant correlators 



Cd- 



CA 



4! 



F A 



CA 



CA 



1 



Computing the diagrams gives 
< AAQaQa > 3 = (l - \ 



-id d p [ 8V2c d c A FA 
(2tt)<* \(M 2 A -p 2 )(M 2 -p 2 ) 
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^(M|-p2)2 (M 2 

= -6V2c A c d F A (l Ml M\ MX 

16vr 2 y 6 g A 2 M 2 - M| g M| y 

I 3 ^ (I I Inr M * 1 Ml I Ml(Ml = ^ lnr M ^ 
+ 16vr 2 1^6 + S A 2 + M 2 - M 2 S + (M 2 - M 2 ) 2 g M 2 ^ " 

We can now add this result to the contribution (|49[) proportional to c 2 ,, and the consistency 
condition translates into a simple relation between the coupling constants Cd and c a 

c A F A = V2c d . (89) 

This allows one to obtain a consistent determination of the correlator component < AAQaQa >3 
associated with the scalar mesons, 

n 3 = < AAQaQa > 3 = — ^ d T A 9N ( i - ^ 9 log ¥4) . (90) 

We can compare the above determination of the coupling ca with experiment. Let us 
consider the decay process /i(1285) — ► ao(980) + tt. Using the Lagrangian above (|S5|l the 
expression for the decay width is easily derived, 

r, J /cm (M\ -Ml)* 

Using the value c<2 ~ 32 MeV from ref.jSS], Ma = \[2M p (and physical masses in kine- 
matical factors) one obtains 

T fl ^ ao7T ~ 7.7 MeV. (92) 

This result is in rather good agreement with the experimental one r^^ ao7r = 8.6 ± 1.7 
MeV from ref . [2] . 

5.6 Contributions from £ wp7r in the vector v/s tensor formalism and LT4 

We discuss here contributions from C wpn ea. ()69|) which are likely to be important and 
were also not accounted for in the minimal resonance Lagrangian. One faces a technical 
difficulty that the Lagrangian ea. (|69|) involves the Levi-Civitta e tensor and this makes 
continuation to d dimensional space ambiguous. One would therefore like to impose a 
stronger version of the consistency condition that the integral be finite in four dimensions. 
This condition is rather simple to enforce if one uses the vector formalism. Let us first 
briefly review what would be the result in this case and then we will show how to proceed 
in the case of the tensor formalism. An extensive study of resonance Lagrangian terms in 
the vector formalism can be found in ref. [IS] • In this formulation, the terms including the 
ujpir type couplings read (we use the same notation as [33]), 

1. „ - , J_ 

2^2" 



+ J^_ e ^ {{u ^ Vv }U a p) + 9 f ^<K, V u }V a p) , (93) 
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where 

Vfj, u = V/jVu - . (94) 

In this description the consideration of short distance conditions imposes that chiral 0(p 4 ) 
terms be added [HH] but these do not play any role for our calculation. Computing the 
correlator < AAQyQy >3 from this Lagrangian, one finds 

f d 4 p (P 2 9a{3 ~ PaPp) \p 2 (9i + 2g\ fv) ~ M v gX] 2 
< AaApQvQv > 3 = 8 / — ^ =r= (95) 



(2vr) 4 p 2 {M v -p 

The integral in eq. (|95|) can be defined in four dimensions without difficulty, it suffices that 
the following relation holds among the coupling constants 

g\ + 2gXf V = • (96) 
The result for the correlator then reads, 

< AAQyQy >r 7r =^(3g 2 v ), ~9V = ^ L - (97) 

The value of gy was estimated in ref.[7] from the ui — > ttj decay rate, \gy \ = 0.91 ± 0.03. 
Under the assumption of nonet symmetry (which is rather well satisfied experimentally 
for the vectors), one also has 

< AAQyQy >2 pn = \ < AAQyQy >^ p7T . (98) 

From eas. (j97|) . 1)98(1 one can deduce an estimate, which is identical to the one performed 
in ref.[7j, for the combination of K^s which is involved in the radiative corrections to the 
process 7r° — > 27. 

Let us now reconsider the problem using the antisymmetric tensor formulation, i.e. 
starting from the Lagrangian of ea.(|69|). One has three kinds of diagrams to compute 
which are shown below 

F v (ft Fy 

(1, 
dj 
Fy 

Evaluating these various diagrams we obtain < A a ApQyQy >3 in the form of the following 
four-dimensional integral 

^ A a n n ^ f d4p I ^PP f MyFyjds + di) r- ^\ 

< A a ApQyQy > 3 = J j-y 4 I y M ,_ p2 - V2(cr + or) j (99) 





32(p 2 g a p - p a p, 3 ) ( 2M v F v d 3 _ r- _ 
M v (M v -p 2 ) [m v -p> VAC5 C6) 



2 
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The integral in ea. (|99|) does not converge in four dimensions unless all the couplings Cj 
and di vanish identically. In order to remove the divergence, in this case, we introduce a 
more general set of resonance Lagrangian terms, which include four building blocks 

C 4bb = {x^sT 1 * + x 2 g^g^ + x 3 g^g^)({f+\ o , U ^}{f^ 2 , U ^}) 

+(z l9 ^g^ + z 2 g^g^ +Z39^ 3 9^ 2 )({Vi ,u tll }{V Xl , 2 ,u tl3 }) . (100) 
These terms generate the following additional contribution to the correlator 



{A a A l3 Q v Q v ) 3 = J 2ip 9aP - p2{M 2_ p2) + p 2 {M 2_ p 




2\2 



~ lttf Ha? ~ PaPp) ( ^ - p2( ^ p2) + p2( M|-V)0 / ' (101) 

In this expression we have used the notation 

£123 = xi + 2x 2 + x 3 , x 13 = x 1 +x 3l (102) 

and similarly with tc%. 

The convergence constraint in four dimension can now be satisfied: it imposes C5 + C7 = 
and determines the values of the six combinations of parameters which appear in eq. (jl01|) 
in terms of d 3 , d 3 + d± and C5 — cq. Collecting all the pieces, one finds that all the couplings 
drop out except d 3 and the correlator < A a ApQvQv > 3 is then expressed as a convergent 
integral, 

< A a A pQvQv > 3 = 128M^§ / ^ ^-^ 8 . (103) 

This expression is exactly analogous to the one obtained upon using the vector formalism 
and we write the final result in the same form, 

F 2 

n 4 = < AAQ V Q V > 3 = — ^(3g 2 AT ), g AT = 4F v d 3 /F . (104) 

l07T z 

The value of d 3 may be determined by enforcing asymptotic conditions concerning the 
< VVP > correlator. Not all conditions can be satisfied at the same time if one uses a 
minimal resonance model [4*21 1431 !"""""] . If one chooses to enforce the overall scaling behaviour 
one obtains |4*T] 

scaling = J% _ ^cM^ ( } 

Numerically, this gives gAT = —0.58 which is somewhat smaller in magnitude than gy . 
Alternatively, one may choose instead to impose that the < VVP > correlator satisfies 
the VMD property, in which case one gets 

rVMD _ N c My 

d3 ~ 6Att 2 F 2 • (1U6) 

If one uses this latter determination, one gets c/at — —0.93, i.e. the result from the 
antisymmetric tensor formalism becomes essentially identical to the one from the vector 
formalism. 
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6 Numerical results 

6.1 Results for K\, K 6 



10 3 K{ 


WK r 2 


10 a Kl 




lQ 6 K r b 


10 3 Kl 




-4.44 


-1.73 


4.44 


-3.46 


13.31 


5.19 


p, ai 


-6.8 


-2.7 


6.8 


-5.5 


20.3 


8.2 


P, ax (BU ) 


1.97 


1.97 


-1.97 


3.93 


-1.97 


-1.97 


copir 


-0.35 


0.35 


0.35 


0.69 


0.35 


-0.35 


scalars 


0.4 


-0.4 


-0.4 


-0.7 


-0.4 


0.4 


scalars (BU ) 


0.11 


0.11 


-0.11 


0.21 


-0.11 


-0.11 


6i(1235) 


-2.71 


0.69 


2.71 


1.38 


11.59 


2.77 


Total 



Table 1: Numerical results obtained for Kl(p) with p, = 0.77 GeV from the present work. 
The contributions associated with the 6i(1235) are discussed in appendix B. The results 
from ref.\2y\j are also shown for comparison. 

In the numerical applications we will simply use M v = 0.77 GeV and F v = V2F n 
(KSRF relation [4*7]). which corresponds to z = 2 (the same values were used by BU). 
We will also use F = F w = 92.4 MeV. The results for K{(p), K%(/i) with p = M v 
are collected in Table 1. We show separately the contributions due to the scalars (terms 
proportional to Cd, with Cd = 32 MeV) and the contributions from the ujpn interactions 
(terms proportional to gAT-, with \§at\ = 0.93). We find that this contribution can be 
significant for some of the couplings like K\ or K\ . The contributions related to the scalars 
are rather small. In that case, our results agree on the magnitude but have a different 
sign from those of BU because of the extra contribution needed to ensure consistency. 
The largest contributions are those originating from the aipir sector. In that case our 
numerical results are in qualitative agreement with those of BU. We have also performed 
an estimate of the contributions associated with the b\ meson (see appendix B) which turn 
out to be rather small. 

6.2 An update on the corrections to Dashen's theorem 

A reliable evaluation of the electromagnetic part of the K + — K° mass difference is im- 
portant for the determination of m u — rrid- At leading order in the chiral expansion it is 
given by Dashen's theorem 26 



It was pointed out that the corrections of order e 2 m s could be rather largeJlS] such that the 
right hand side could be modified by as much as a factor of two. Using Urech's formalism 
one can compute explicitly these corrections which have the following form [3] 



(AMI)em = (AM^m + 0(e 2 m q ) . 



(107) 



(AM|) EM - {^MDem = e 2 M 2 K (A 1 +A 2 + A 3 ) + 0(e 2 M 2 ) 



(108) 
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with 



Ax 
A 3 



167T' 



■[(3 + 2Z)log 



Ml 



4]-16Z27 5 , Z=^ 



8(K[ + «Ti) 



(109) 



As discussed in (I) (and also in (301) the Urech couplings K{ , K\\ taken separately 
suffer from short distance ambiguities and depend on the gauge. The sum, however, is 
perfectly well defined and can be estimated from sum rules related to the correlators 
< y3y3 _ y8y8 _ ^43^43 + ^8^48 > and < y^p > gi v j n g the following result [3T] 



K[ + K 
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(110) 



(111) 



The quantity Za — Zy is generated from flavour symmetry breaking at first order in 
the vector multiplet as well as the axial-vector multiplet and the computation takes into 
account Weinberg-type sum rule relations for the splittings of the masses as well as the 
coupling constants. The last entry in the expression (|110[) . is a contribution which 
should restore the scale dependence proportional to C/Fq . Since this effect is subleading in 
N c , this term will be dropped here for consistency with previous approximations. Making 
use of the result eq. (|110j) and of the numerical values for K£ + Kq shown in table 1 we 
find the following numerical values for the three corrective terms (which we give in units 
of the physical ir + — 7r° mass difference, AM 2 = 1261.16 MeV 2 ) 



(AMl) EM - {AMDem = AMI [0-59 - 0.34 + 1.25] 



(112) 



The new result with respect to (I) concerns K§ + Kq and we find that this contribution 
goes in the sense of reducing the size of the correction. This is in qualitative agreement 
with the result of BU. However, the effect of the extra terms from C ai(m included here 
is to reduce the magnitude of K§ + by approximately a factor of two as compared to 
BU. There is a clear indication for a sizable correction to Dashen's theorem rather similar 
to that found in refs. |49| I3()j. The origin of the different result obtain in ref.[2§| can be 
traced, essentially, to the fact that their result for K\q + Kh is ten times smaller than ours. 
This is due to two reasons 1) flavour symmetry breaking in the vector and the axial- vector 
multiplets was ignored in [21] such that K\q was just set to zero and 2) the effect of the 
couplings from C aipn which they did not take into account is to strongly enhance the value 
of K n . 
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6.3 Uncertainties 

One should keep in mind that there are uncertainties which are difficult to estimate quan- 
titatively. These have two main sources. The first one is that we have used, at several 
places, large N c approximations a) in the treatment of the resonances, which are taken 
to be infinitely narrow and b) in the calculation of QCD n-point correlators in which 
only tree graphs were considered. The second source of uncertainty is that we have re- 
tained the contributions from the lightest resonances only. The usual experience with such 
approximations is that reasonable order of magnitudes should be obtained. 

7 Conclusions 

In this paper, we have considered a set of sum rules for the electromagnetic chiral pa- 
rameters which involve QCD 4-point correlators in connection with the construction of a 
chiral Lagrangian with resonances. For our purposes, the resonance Lagrangian provides 
rational-type approximations for the QCD correlators obeying the chiral Ward identities. 

We have shown that the resonance Lagrangian of order p A was not sufficient to ensure 
all the convergence conditions required by the sum rules. We have then shown that 
introducing a set of terms of order p e involving aip7r-type couplings renders it possible 
to satisfy these conditions in the sense of finiteness as a d — > 4 limit. The new coupling 
constants which appear were shown to obey non linear equations, we have argued that 
physical requirements select a unique solution. 

Our investigation, however, was essentially limited to a set of terms with three building 
blocks. A more general investigation should be performed in the future. A more complete 
determination of the resonance Lagrangian would lead to much better estimates of the 
chiral coupling constants at this order [SOj than those available at present. This would 
improve significantly the effectiveness of the chiral calculations at order p 6 . There is also 
a limitation in the set of resonances which were included. In this regard, the role of tensor 
mesons should perhaps be investigated |51). 

We have considered an application to the problem of the chiral corrections to Dashen's 
low-energy theorem. We have found that the combination + K£ as determined here 
reduces the size of the correction which still remains rather large. One motivation for 
deriving estimates for each of the parameters Ki is to improve the accuracy of radiative 
correction evaluations. In this context, further work is still needed because another set of 
parameters, JQ, is involved in the interesting case of the semi-leptonic decays [Sj. Finally, 
a formal similarity can be noted between the sum rules considered here and those shown 
recently to hold for chiral parameters in the weak non-leptonic sector [52]. 

A Vertices 

We list below the set of vertices which are needed in sec. 2 and those needed sec. 3 which 
are associated with vector and axial-vector resonances. The vertices represent derivatives 
of the action (times i) with respect to the various sources (with each source multiplied by 
i). It is not difficult to derive, in an analogous way, all the other relevant vertices. 
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A.l Vertices with pions 

V = -Fd x n c (z)F x (z) 



V = iFf cat 7r t (x)F a (x)5(x - z) 



V = iFf bct 7r t (y)F f3 (y)5(y - z) 



V = iFf bat Tr t (x)g a ^5{x - y) 



V = -f stc 8^ s (z)n\z)F^z) 



V = -f sta d a TT s (x)it\x) 



v a a {x) X» q c A (z) 



a b 3 (y) + 



O q c v {z) 



v^x)x+a b fi (y) 



V = ig a Pf ais f b i s ir i (x)iri(x)5(x-y) a a a {x) 4 



4(y) 



V = -ig a ^f ais f b ^ s 7r i (x)^{x)5(x - y) v a a {x) XX v b {y) 
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A. 2 Resonance vertices with sources: 

These vertices are the same for vector or axial-vector mesons: 



V = -f^R^{x)d x R^{x) J v a a (x) 



V = if ail f h ^ l R i 0lll (x)Rf{x)5(x - y) v%(x) XX v h p (y) 



A. 3 Vertices proportional to F A 



V = —FaA c ^ v (z) d fl F u (z) 



V = iF A f bc AUy)F\y)5(y-z) 



V = iF A f iad A i aX (x)F x {x)5{x) 



V = ^f^A* lv (0)& l F v (0)^(0) 



• <f A ( z ) 



a a a (x) +o^(0) 



4(0) ol 



Note that this vertex is antisymmetric in the Lorentz indices a(3 and will therefore not 
contribute here. 



36 



A. 4 Vertices proportional to F v 



V = F v V^ v {z)d^F v {z) 



O q v (z) 



V = -iF v f mc V^ x (x)F A (x)5{x - z) 



v a a {x) XO q c v (z) 



V = -iF v f iac V^ x (x)F x (x)5(x - z) 



V = 2iG v f act V^(x)F u (x)5(x - z) 



a a a (x) +• q c A (z) 



V = ^ff StC V^F U TT t 



A. 5 Vertices proportional to G\ 



a%{x) +• q c A {z) 



V = =™x-f'»*V°(x)d»ir t ( 



V = ^§^f cst V^(z)d^TT t {z)F x ( 



V = =^f stu j dvV"(v)d^ t {v)d v -K u (v) 



a%{x) +1 



q c A (z) O 



There are other potentially relevant vertices which are antisymmetric in a(3 and will 
thus not contribute for our purposes. 
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B 6i(1235) contributions 



We present here an estimate of the contributions associated with the b\ (1235) resonance. 
The physics of the b\ meson is not as constrained as that of the a\ meson by chiral sym- 
metry, so in this section we will occasionally make some educated guesses. The order of 
magnitude of the b\ effect should nevertheless be properly accounted for. The contribu- 
tions are to the same correlator components as those from the uopix couplings. In principle, 
the convergence constraints should be applied to the sum. We will applyhowever, these 
constraints separately to the b\ for the purpose of obtaining an order of magnitude esti- 
mate. 

We can write down a set of vertices in exact analogy with the aipn vertices 



(GfgWg™* + G B g^g^ + G B g™*g* l >»){ {V M ^ lA , B*}u 



"M3/ 



+ ( H B g twi g M3 + H B g wm g m»i + HggWSg^^dV^x, B A }V 



(113) 



The only difference with the Lagrangian of eq.(|67|) is that the commutator is replaced by an 
anticommutator on account of the fact that the b\ meson is odd under charge conjugation. 
As a consequence, the contribution to the correlator {AAQyQv)'i can be deduced simply 
from the results of sec 15 .41 



{AAQvQ 



biUJTT 



+ 



-3F 2 
16vr 2 

7M 2 - 



(Gf 



G B + 



2G B \ 



log 



M 2 



r4 1 A£S 



M%\og 



A 2 (M 



My) 2 



Ml 



6(M 2 



Ml) 



+ 4 



(Gl) 2 My 



Ml 



l0g A2- + 3 



(114) 



Imposing a consistency condition on eq. (|114j) enables one to express (AAQvQv)z 
terms of the combination 



b\WK 



111 



G 



Gf + G B + 2Gf 



and one obtains 



(115) 
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We must next try to determine the relevant combination of couplings G B from experiment. 
The widths of the decays bf — ► jtt + and bf — ► luit + are known 



+771" 1 



227 ± 59 KeV, T h 



142 ± 9 MeV 



(117) 



and their expressions in terms of the resonance Lagrangian coupling constants read 
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Here, we have introduced the notation 

H B = H B + H B + G B + 2G B (119) 

Using these two pieces of information enables one to determine the combination of cou- 
plings that one needs, but there are several solutions due to the non linearity of the 
equations 

H B = ±0.68 , G B = Gf + G B + 2G B = ±0.24, ±0.67 . (120) 

We can eliminate some of the solutions by using one more piece of experimental information 
concerning the D/S ratio in the decay b\ — > lott 

^ = 0.29 ±0.04 (121) 
Js 

The S and the D-wave components of the decay amplitude b\ — ► coir (defined as in ref . ) 
are easily computed from our Lagrangian 

B 



fs = ~ zm b M^ ( Mb " [ Mv ^ Mb + 2M v)G B + {M 2 B - Ml)R 



f*> = - fi,/^ f ( M b ~ M v) \2M V (M B - M V )G B ± {Ml - My)H B ] (122) 

Using this result and the experimental one on the D/S ratio reduces the set of acceptable 
solutions for the coupling constants to 

{G B ,H B ) = ±(0.24, 0.68) (123) 

which give fo/h = 0.28. 
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